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Abstract 

The interaction between planar quasiparticles in a d x 2_ y 2 superconductor and quan- 
tized vortices associated with a magnetic field perpendicular to the plane is shown 
to induce a pair potential with d xy symmetry, out of phase with d x 2_ y 2 order. A 
microscopic calculation of a process involving quasiparticle scattering by the super- 
current around a vortex and Andreev reflection from its core is presented. Other 
processes also leading to an id xy pair potential are discussed. It is argued that such 
a fully gapped state may be the high field low temperature phase observed by Kris- 
hana, Ong et al in magnetothermal conductivity measurements of superconducting 
single crystal Bi — 2212. 



I. Introduction 

Recent measurements of in-plane thermal conductivity Kt in superconducting single 
crystal Bi — 2212 (1,2) for T ^ T c and H <C H c2 show a rather sharp change in its 
magnetic field dependence at a field H (of order a Tesla or more) which depends on 
temperature T Q (~ 10K) nearly quadratically. For H > H Q and T < T a , K T hardly 
changes with the magnetic field, in sharp contrast to the decrease exhibited below 
H Q . This suggests (1,2) a transition from a d x 2_ y 2 phase to a fully gapped phase at 
H (T ) with an exponentially small density of heat current carrying quasiparticles. 
I show here that in a magnetic field, a d x 2_ y 2 superconductor necessarily develops a 
pair amplitude of id xy symmetry at low enough temperature, leading to a gapped 
phase which might be the basis of the observed Kt behaviour. 

In the superconducting state, which has vortices due to the external magnetic field, 
the gap quasiparticles interact with the circulating supercurrent around each vortex, 
as well as with the order parameter inhomogeneity associated with the vortex core. 
These interaction terms are obtained explicitly (Section II). The first scatters quasi- 
particles, eg. it changes their phase. The second causes Andreev reflection. Due 
to their combined effect, a particle interacting with a vortex comes out as a phase 
shifted hole, i.e. the vortex is a source of out of phase pair potential. This term is 
calculated in Section III, and is shown to have id xy symmetry. Thus, at each vortex, 
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id xy order is necessarily induced. If the vortices are well ordered spatially, there is a 
homogeneous A xy term in the free energy, proportional to the density of vortices or 
the magnetic field B. I argue (Section IV) that this could be the phase suggested by 
the observations of Krishana et al (1,2). Related results and questions, such as the 
temperature and field scale of the transition, effect of vortex lattice periodicity are 
also discussed (Section IV). Recently, Laughlin has proposed (3) that in a magnetic 
field, the pair potential is of the form d x 2_ y 2 + a id xy where a is real, and that 
there is a (first order) transition to such a phase in the (H, T) plane. The proposal 
is based on a mapping of a d x i_ y i + a id xy Hamiltonian to a lattice quantum Hall 
system, and on exploiting quantum Hall current ideas. Here, I describe a detailed 
microscopic mechanism for id xy order, and calculate its size. 



II. Quasiparticle Vortex Interaction 

I show here that the quasiparticle Hamiltonian in the mixed state can be written 
as a sum of two classes of terms; one is that of the uniform superconductor, and 
the other describes the effect of vortices on quasiparticles. Since quasiparticles are 
well defined in the superconducting state, a mean field Hamiltonian is adequate. 
In the homogeneous superconductor, the otherwise free electrons move in a real, 
spatially homogeneous pair potential A with d x 2_ y 2 symmetry, eg. with Aj, = 
A G (cos k x a — cosk y a) for nearest neighbour pairing. In the presence of vortices, 
the pair potential A(r, ri) is inhomogeneous and complex. Its phase changes by 2n 
on going round a vortex, and the magnitude must vanish at each vortex core Rg. 
Separating out the phase part, we can write A(r, ff) as 

A (fff) = A m (f - ff, R-R e )exp[ l -{J2 0(r-R e ) + 9 (ff - R e )} ]. 

1 i 

Here R is the centre of mass coordinate if + fl)/2. Making a gauge transformation 
ip^(r) exp{ | J2e @(f — R^} =>- ip^(f), the quasiparticle Hamiltonian becomes 

H = / drip„{r) (p+ mv s (r)) 2 i{j a (r) (1) 

+ J dfdff[A m (f- ff, {R - Re} )^{r)^{fi) + h.c] (1) (2) 

where 

v s (r) = (2m)" 1 ]T {hy9 (f - R t ) - 2eA(f- R e )/c} 

i 

is the in-plane gauge invariant superfluid velocity. In the London limit, A m can 
be considered spatially uniform [R independent) with d x 2_ y 2 symmetry, except for 
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order parameter deficits 5A m around each vortex. (These have to be determined 
self consistently) 

The quasiparticle Hamiltonian Eq.(l) can be written as the sum of an unperturbed 
term H Q , and the remainder describing quasiparticle vortex interaction. In the 
momentum representation, and the Nambu formalism, one has 

H = H + (H e + H m + H KE ) (4a) 

where 

H o = J2 a l (efc r 3 + A fc T! )a k (46) 

k 

H e = ^a\Tik.v s (q)a k _ q (4c) 

k,q 

H m = -A fkgialn a k - q ) (Ad) 

k,q 

and Hke is the superfluid kinetic energy. Eq.(4b) describes a d x 2_ y 2 superconductor. 
He is the quasiparticle superfluid velocity interaction and H m is the inhomogeneous 
pair potential due to the vortex core, with — A /^^. being the Fourier transform of 

the deficit SA m (r — rf, {R — R^} ). The superfluid velocity v s (q) is assumed to have 
the standard Ginzburg Landau form 

away from the vortex cores. For wavevector transfers q <^ A^ 1 , this has the un- 
screened form independent of A 2 . We now use the mean field Hamiltonian Eq.(4) to 
show that id xy order is induced in a magnetic field. 



III. id xy Order 

The form Eq.(4) is natural for looking into the question of the pair potential in the 
presence of a magnetic field. In its absence, HQ,H m and H KE vanish, and A k in 
Eq.(4b) is just the d x i_ y i uniform value. We imagine H Q and H m being turned on, 
and ask if any out of phase order is induced, i.e., in the notation above, whether 
Afc = (a k T2dk) 7^ 0. We are also interested in the k dependence of \ k . We shall 
calculate A& with Hq and H m as perturbations. This enables us to focus directly on 
processes leading to a particular order parameter symmetry. 

The first nonvanishing contribution to A& comes from a term linear in Hq as well 
as in H m ; there are no contributions to second order, either in Hg or H m . The 
process is diagrammatically represented in Fig. la. It describes the combined effect 
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of quasiparticle scattering from state k to (k — q) by the supercurrent due to the 
vortex, and Andreev reflection from (k — q) to a hole of momentum-A; caused by the 
inhomogeneous order parameter magnitude associated with the vortex core. 

The pair amplitude A^ can be calculated from Fig. (la), and Eq.(4). At T = 0, it is 

{a\r 2 a k ) = £(A / fc? ){ {K 2 /2mA) (k x iq.e z )<T 2 } [(e fc _,/ ' E k E k „ q {E k + E k _ q )\. (6) 
i 

In Eq.(6), the first bracket is the pair potential due to the vortex core; approximately, 
fkq = fq{cos(k x + q x )a — cos(k y + q y )a} where J2 q fq = 1- The second factor, in 
curly brackets, is the quasiparticle supercurrent coupling k.v s (q) for q ^> A -1 (the 
regime of interest). The last term arises from the intermediate state sum, with 
E k = |(e| + A!) 1 / 2 ! • On expanding the integrand as a power series in q, the leading 
term is 

\ k = a (n v /n) (A /E k ) (k x k y ) (7) 

Here, a is a constant of order unity, and n v is the vortex density. 

We notice that A^ has d xy symmetry. It is a result of the coupling between the 
internal and centre of mass states of the pair near a vortex. In the Fourier represen- 
tation used, the former is described by the momentum k, and the latter by q. The 
quasiparticle supercurrent coupling has a structure (k x q).e z that affects the angu- 
lar state of a pair, through the constituent single particle dispersion. Algebraically, 
(k x q).e z {^ k _q) — k x k y (g 2 — g 2 ). This term is clearly nonlocal, ie. it arises from 
higher powers of q or of gradient/derivative terms. Thus, a fully local, semiclassical 
theory, working in terms of a spatially slowly varying superfluid velocity v s (r) and a 
local (diagonal) quasiparticle momentum or Fermi surface (4), will miss this effect, 
while it may be appropriate for quasiparticle density of states etc.. The id xy order 
can also be induced by higher order terms involving only Hg or the quasiparticle- 
supercurrent interaction. An example is shown in Fig. (lb). Other contributions, 
involving induced s wave order, are also possible. 

We have focussed so far on spatially uniform single site terms. The id xy order induced 
around each vortex is spatially nonuniform. The leading nonvanishing contribution 
is of first order in H e , and is shown in Fig.lc. It describes a d x 2_ y 2 pair becoming a 
d xy pair and simultaneously acquiring a momentum q, near a vortex. For q ^> A -1 , 
the contribution from this diagram is nonvanishing at small q. However, because of 
London screening, the q A -1 limit vanishes. A related process involving d xy order 
near a spin orbit impurity potential has been recently discussed by Balatsky (5). 

The perturbative approach described above raises the obvious question of the ex- 
pansion parameter. On general coupling constant and phase space grounds this can 
be argued to be (l//cp£) for both H 9 and H m , where £ is the pair coherence length. 
This is about (1/5) or (1/6). Thus, perturbation theory is expected to converge. 
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The id xy order Xk leads to a gap iA xy with xy symmetry, in two ways. A process 
for Afc corresponding to Fig. la has, associated with it, an anomalous self energy 
Fig. Id, which leads to an xy symmetry gap function. Thus even if there is no 
interaction in the xy particle particle channel, a A xy is induced. Secondly, there may 
actually be such a potential (attractive or repulsive), say V xy . Then in the mean field 
approximation, a gap A xy (k) ~ J2kr Vxy (kk') X^ is induced. The question of self 
consistency can be fully addressed only if the microscopic mechanism of pairing (say 
in the d x i_ y i channel) is known. In a BCS theory, if appropriate pseudopotentials 
V x 2_ y 2 and V xy are assumed, the relevant Gor'kov equations (in the presence of 
vortices) or the quasiparticle problem (described by Bogoliuov-de Gennes equations) 
need to be solved self consistently. We have not done this; the lack of self consistency 
does not affect either the existence or the approximate size of the T = id xy pair 
amplitude. 

IV. The Gapped Phase 

We have shown above that an id xy pair amplitude is inevitable for a d x 2_ y 2 superconductor 
in a magnetic field, and that there is a consequent id xy gap. Thus, the ground state 
of the system in a magnetic field is necessarily a ful ly gapped superconductor in 
which the gap parameter = ^ A 2 x2 _ y 2(k) + Al y (k) is nonzero at all points k on 
the Fermi surface. This has obvious thermodynamic and transport consequences. In 
particular, since the number of quasiparticles is exponentially small at very low tem- 
peratures, their contribution to thermal conductivity is negligible. We thus expect, 
at low enough temperature and sufficient vortex density, that there will be a tran- 
sition to a gapped phase. We do not have a complete theory of this transition, but 
discuss possibilities below, following a brief analysis of experimental results (1,2). 

As mentioned earlier, Krishana et al (1,2) find a transition at T Q from a K T which 
decreases with increasing if to a field independent value at H ; approximately, 
T Q oc \[TT . It is argued that this implies a thermodynamic transition to a gapped 
or coherent phase from an ungapped, incoherent phase. The transition is most 
likely continuous, since for a discontinuous change one expects a jump in K T at the 
transition point. 

The physical picture of the transition is that around each vortex an inhomoge- 
neous d xy order develops. If the vortices are not regularly arranged, and if the d xy 
amplitudes are small as well as patchy or disconnected, quasiparticles see an inho- 
mogeneous medium which scatters them, and a pseudogap develops near the nodes; 
this pseudogap deepens as the vortex density increases so that the electronic ther- 
mal conductivity decreases with increasing vortex density. At some critical vortex 
density dependent on temperature, the id xy order parameters overlap, or the vor- 
tices order spatially, and a nonzero gap develops for the lowest lying quasiparticle 
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excitations. This is the new phase. 

If the d xy pair amplitude, (labelled m^ y ) is the expected order parameter, an order 
parameter functional can be obtained using the standard auxiliary field method. 
We have calculated low order terms at T = in such an expansion. There is a 
linear term, due to the fact that averages (aj^Ofc) with k x k y symmetry, discussed in 
Section III, are nonzero. There are quadratic terms arising from both V xy symmetry 
interactions between electrons (pairs) as well as a nonzero xy pair susceptibility. 
There is a cubic kinetic energy term (3). Because of the linear term, to* is nonzero; 
unsurprisingly, it has the value (a'^a^xy for which some terms have been indicated 
in Section III. We have not carried through the calculation for T ^ 0, and so are 
unable to find the temperature below which m k xy becomes nonzero. The linear in to 
term persists however at T ^ 0. 

The gap in the excitation spectrum at T = is due to both anomalous self energy 
(Fig. Id) and the mean xy pair potential if there is a V xy interaction. The gap due to 
the former is approximately A (n v /n) and due to the latter, 4 ((V^l/e^) A (n v /n). 
Here A G is the d xy gap magnitude and n is the carrier density. The actual value 
thus depends on V xy which is not known. Assuming V xy /ep ~ 1 (a large value) and 
A G ~ 300i^, for a field of 5T, the minimum gap A xy is about 10K. This comparable 
to the temperature 20K at which the gapped phase transition occurs for 5T, though 
smaller. Also, equating T D to A xy , we find that T Q oc H; the experimental points are 
closer to y/~H. 

This suggests that the identification of the observed transition with the development 
of a uniform part of d xy pair amplitude may not be correct. One curious feature of 
this identification is that even when the vortices are distributed randomly the pair 
amplitude is uniform (there is a diagonal in k term) It is the sum of N v independent 
identical single vortex terms irrespective of their arrangement. However, the d xy 
order developed due to quasiparticle vortex interaction is spatially inhomogeneous, 
on the scale of the screening length A, as discussed in Section III. If the vortices 
order spatially, there are coherent terms linear in H 9 , with wavevectors q equal to 
Bragg vectors of the reciprocal vortex lattice. (These are q = G components of the 
term shown in Fig.lc). They have the right energy scale. The relation between 
the magnetic field and the related energy scale is T c oc G or T c oc \/~H, close to 
that observed. The gapped nature of this vortex lattice phase, the preferred lattice 
structure, and the nature of the transition to this lattice on cooling, are being 
investigated (6). 
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Figure Captions 

Diagrams describing processes leading to id xy pair amplitude and related anomalous 
self energy. 

(a) An interference term for id xy order involving a quasiparti- 
cle (straight line) being scattered by supercurrent (wavy line) 
around a vortex (marked by a cross), and by order parameter 
inhomogeneity around the vortex core (dotted line). 

(b) A higher order term for id xy order, solely from quasiparticle 
supercurrent coupling H Q . 

(c) A spatially inhomogeneous id xy pair amplitude term, first order 
in H e . 

(d) Anomalous self energy connected with the process of Fig. (la). 
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